Local Electronic Correlation at the Two-Particle Level 
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Electronic correlated systems are often well described by dynamical mean field theory (DMFT). 
While DMFT is an exact self-consistent theory in the limit of infinite coordination number or 
dimensions, when applying it as an approximation to finite-dimensional systems, its self-consistency 
is guaranteed at the one-particle level only, i.e., for the local Green's function and self-energy. 
However, extremely valuable information is also enclosed in the local two-particle Green's functions 
and vertices. In fact, even at the DMFT level, mastering the calculation of two-particle quantities 
is crucial to compute momentum-dependent response functions that can be compared directly with 
experiments. Moreover, the knowledge of the local two-particle Green's function represents the main 
ingredient for including non-local spatial correlations at all length scales by means of diagrammatic 
extensions of DMFT, such as the dynamical vertex approximation or the dual fermion approach. 
As hitherto the investigation of local two-particle properties has been merely sporadic, we present 
here a systematic analysis of the local reducible and irreducible two-particle vertex functions for the 
Hubbard model at the DMFT level in the context of an unified diagrammatic formalism, providing 
an interpretation of the observed frequency structures in terms of perturbation theory results and 
of the mapping onto the attractive Hubbard model. 

PACS numbers: 71.27.+a, 71.10.Fd 



I. INTRODUCTION 



Electronic correlations are responsible for some of the 
most fascinating phenomena occurring in condensed- 
matter physics, such as the colossal magnctorcsis- 
tance of the manganitesi, the Mott-Hubbard metal- 
insulator transition (MIT)2 in the vanadates^, the high- 
temperature superconductivity of cuprates^ and (pos- 
sibly) of iron-pnictides^, and even for the appearance 
of quantum critical points in particular heavy fermion 
compounds^. While the exact treatment of electronic 
correlation is an impossible task in real materials as 
well as in model systems (e.g., the Hubbard model 7 -), 
a major step forward has been obtained since the 
early Nineties with the dynamical mean field theory 
(DMFT^. DMFT represents the quantum extension 
of classical mean field theory, and, hence, can be rigor- 
ously derived as the exact solution of a quantum many 
body Hamiltonian (such as the Hubbard Hamiltonian) 
in the limit of infinite coordination number or dimen- 
sions (d — > oo)§.. While the average over infinite spa- 
tial dimensions implies neglecting non-local spatial cor- 
relations, DMFT provides for a very accurate treatment 
of local quantum (dynamical) fluctuations. In the case 
of localized electrons, these fluctuations play a pivotal 
role, as they can drive, e.g., the MIT in several com- 
pounds. The more convincing proof of the accuracy of 
DMFT and/or its combination with ab-initio methods 
(LDA+DMFT)^, however, comes from its impressive 
success in treating and explaining some of the most chal- 
lenging problems in condensed matter physics. We recall, 
among the most successful applications of DMFT, the de- 
scription of the S phase of Pv±ii, the MIT in V 2 0^, the 
correlation effect in Fc and Nii£, and the volume collapse 
in Ce^. Using DMFT has become almost standard for 



treating electronic correlated systems in the last decade 
and, as an example, recent DMFT calculations have been 
able to explain the appearance of kinks in the spectral 
functions and the specific heat of particular vanadates (as 
SrVOy^, LiV^O^i^), the spin-polaron peak structures 
in photo-emission^ 7 - and optical spectroscopic data^ of 
strongly coupled antifcrromagncts, such as V2O3 and 
LaSrMn04, the anomalies of the optical sum rule behav- 
ior in the high-temperature superconducting cupratesi^, 
as well as in V^O^, and some of the spectral^ and 
magnetic properties^ of iron-based superconductors. 

However, by looking in more detail at the exist- 
ing DMFT or LDA+DMFT studies, one clearly sees 
some limitations. For instance, the theoretical cal- 
culations and the comparison with experiments are 
mostly performed for one-particle quantities only, such 
as momentum-integrated or momentum-resolved spec- 
tral functions. The analysis of the two-particle quan- 
tities is usually restricted to the easiest cases of optical 
and thermal conductivities, for which -in DMFT- one 
can safely consider— the "joint-density of state" ("bub- 
ble") term only, i.e., in other words, it essentially re- 
mains at the one-particle level. Calculations of spec- 
tral properties at the "actual" two-particle level are - 
with few exceptions^—- done for local susceptibilities 
only, which can be also approximated with those of the 
sclf-consistcntly determined Anderson impurity model 
(AIM) associated with DMFT. The reason for this is - 
evidently- the higher difficulty and heavier workload of 
performing calculations at the two-particle level. In fact, 
as it was already reported in Rcf. |f| the calculations 
of a (spin, charge, particle-particle) momentum- and 
frequency-dependent susceptibility x(q,^) at the DMFT 
level requires the determination of the irreducible local 
vertices (r r , with r = d, m, s, t ) of the AIM in the cor- 
responding particle- hole (density/ magnetic) or particle- 
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particle (singlet/ triplet) channel, which serve as an input 
for the related Bethe-Salpeter equation. 

The importance of determining the properties of re- 
ducible and irreducible two-particle quantities for the 
AIM however, goes well beyond the calculation of the 
momentum- and frequency-dependent response functions 
x(q, cj), needed for the comparison of DMFT with other 
spectroscopic experiments than photo-emission or optics. 
Indeed, reducible and irreducible vertices of the AIM 
are the basic ingredients of two important diagrammatic 
extensions of DMFT, such as the dynamical vertex ap- 
proximation (DrA) 26 i 27 and the dual fermion (DF)2& ap- 
proach. In fact, both methods aim at the inclusion of 
spatial correlations beyond DMFT at all length scales, 
starting from a two-particle level (local) input of an as- 
sociated AIM. Leaving aside the theoretical and numer- 
ical challenges of performing the calculation for the lo- 
cal two-particle vertices^, a thorough analysis of their 
general properties and of the physical interpretation has 
been lacking in the literature hitherto. The full frequency 
dependence of two-particle quantities indeed has been 
shown or discussed only in selected cases and for very 
specific problems (e.g., Refs. H^, [2(| andl3ll). The main 
scope of the present paper, hence, is to fill this gap. 
We provide a systematic DMFT study of the two-particle 
reducible and irreducible local vertices of the Hubbard 
model within a unified derivation and formalism. In this 
framework, the interpretation of the main structure of 
the vertices in (Matsubara) frequency space will be made 
easier by the comparison with perturbation theory re- 
sults and by the mapping onto the attractive Hubbard 
model. Both the formal derivation and the physical in- 
terpretations are potentially of high impact for future 
developments of many-body theoretical schemes and for 
calculations relying on an increased understanding and 
inclusion of correlations beyond the one-particle level. 

The scheme of the paper is the following: In Sec. II, we 
first introduce the formal and diagrammatic definitions 
for treating the reducible and irreducible two-particle 
Green's and vertex functions of the AIM associated with 
the self-consistent solution of DMFT. At the end of this 
section, we also mention the general symmetry proper- 
ties that are expected for such vertices (while their formal 
derivation is explicitly given in the Appendix). In Sec. 
Ill, our DMFT results for the reducible and irreducible 
local vertex functions are presented together with their 
interpretation in terms of the corresponding pcrturba- 
tive results. Furthermore, we analyze the effect of differ- 
ent approximations at the two-particle level on selected 
physical quantities. Subsequently in Sec. IV, the map- 
ping onto the attractive Hubbard model is exploited to 
gain further insight into the main structures of the two 
particle vertex functions. Finally Sec. V is devoted to 
summarizing our theoretical and numerical results and 
conclusions. 



II. TWO-PARTICLE DIAGRAMS: FORMALISM 
AND GENERAL PROPERTIES 

Starting point for our analysis is a rigorous and co- 
herent definition of the relevant one and two-particle 
quantities and of their general properties, which we will 
use throughout the present paper. While part of the 
derivations reported in this section (and in the corre- 
sponding appendixes) is already know n 33 ' 34 , to the best 
of our knowledge, a systematic and unified discussion of 
the two-particle properties has been reported only par- 
tially or implicitly in the standard literature of quantum 
field theory of many particle systems. Hence, the ex- 
plicit derivation of local vertex definitions and properties 
is helpful for an easier reading of the following sections, 
where our numerical and analytical results are presented. 

As mentioned in the Introduction, we consider one of 
the most fundamental models for electronic correlations, 
the Hubbard model on a simple cubic lattice 

^Hubbard = ~t ^ cjcAv + U ^ (1) 

Here t denotes the hopping amplitude between nearest 
neighbors, U is the on-site Coulomb interaction, and 
& ia [cia) creates (annihilates) an electron with spin a on 
site i; hi„ = c ia Ci a . In the following, consistently with 
previous DMFT and DrA paper a 26 ' 35 , we will express 
all energies in units of D = 2\/6t, which ensures that 
the standard deviation of the non-interacting density of 
states (DOS) is kept fixed to 0.5^. As we are mainly in- 
terested in purely local quantities at the DMFT level, for 
indicating these we will omit the site-index i in the fol- 
lowing. Specifically, as in the DMFT-limit of infinite co- 
ordination number, the Hubbard-model can be mapped 
onto an effective (self-consistently determined) AIM. We 
will use the latter to calculate analytically or numerically 
the local observables 

H = Haim = eia\ a ai a +y^ j Vj (c\ ai a +a\ a c a )+U hfhi , 

la la 

(2) 

where a\ a {ai a ) creates (annihilates) an electron with spin 
a at the bath- level of energy et, cJ(c CT ) creates (annihi- 
lates) an electron at the impurity site (n a = c\c a ) 1 Vi de- 
scribes the hybridization between the bath and the impu- 
rity, and U is the on-site repulsion between two electrons 
at the impurity. 



A. Definitions and general properties 

The general definition of the /(-particle Green's func- 
tion G n reads^ 3 - 

G n ,ai...a2 n (n, ■ ■ - ,r 2 „) := (T(c] ri (r 1 ) . . . c CT2ri (r 2n ))), 

(3) 



3 



where an odd/even index always corresponds to an cre- 
ation/annihilation operator $ a /c a . This means that the 
creation and annihilation operators appear in alternat- 
ing order in Eq. and (6) = ^Xx{eT m 6) with Z = 
tr(e~^ n ) denotes the thermal expectation value for the 
observable O. T denotes the time-ordering operator—. 
For n = 1, one obviously recovers Gi i(Ti(T2 (ti,T2) = 
Go-{ti,T2) = G(ti,T2), i.e., the one-particle Green's func- 
tion. In the two-particle case (n = 2), one usually con- 
siders the so-called "generalized susceptibility", defined 
by the following combination of one- and two-particle 
Green's functions 



X<7\<7l (T3CT4 

(Vl,T2,T3,T 4 ) := G 2 ,ITi...IT4 (n, 

— Gi !tTia2 (n, r 2 )Gi iCT3Cr4 (r3, T4). 



(4) 



Without any loss of generality, one can always limit the 
domain for the imaginary times tj to the interval [0,/3] 
(cf. Ref. [IH and Appendix |A"]) . Furthermore, due to the 
time-translational-invariance of the Hamiltonian, one can 
restrict oneself to only three time-arguments T\, t 2 , t 3 in 
the interval [0, /?] (cf. Appendix[Aj i.e., we can set T4 = 



X 



U l (7-2 &30"4 



(n,T 2 ,r 3 ) := G 2iCri ... cr4 (Tl,T 2 ,T 3 ,0)- 

"Gi i<Tl(J2 (ri, r 2 )Gi !(T 3 (T4 ( T 3> 0). 



(5) 



particle (pp) notation, respectively 



A J ph,<7<7 1 



x( vo, {y 1 4- bS)a [y + := 

outgoing electrons incoming electrons 



dTidT 2 dTz Xaa' (ti,T 2 , T 3 )e 



-iwri i(i/+u>)T 2 p -i{v +w)t 3 



(7a) 



outgoing electrons incoming electrons 



/ dry dr 2 dr 3 Xaa' (n , r 2 , r 3 ) e 
Jo 



(7b) 



with z/ and i/ being fcrmionic Matsubara frequencies (i.e., 
v, v' = jj(2n +l),n £ Z) and w being a bosonic Matsub- 
ara frequency (i.e., ui = S(2n), n £ Z). 
The choice of the frequency convention for both cases 
has a clear physical motivation, (i) In the ph-c&se one 
considers the scattering process of a hole with energy — v 
and an electron with energy v + w, i.e. the total energy 
of this process is u. 



One should also recall that, for the SU(2)-symmetric case 
considered here, the spin-indexes a\ . . . 04 are not com- 
pletely independent, as a results of the conservation of 
spin. In fact, among the 2 4 = 16 possible combina- 
tions of spins, only the following 3x2 = 6 remain: (i) 
o\ =0-2=0-3 = cr 4 , with cri=t,|; (ii) (01=02)^(03 = 0-4), 
with ui =ti^; (iii) ("I = °4) ^ (02 = 03), with a\ =t,i. 
This suggests the following definitions 



X<T<T'(n,T2,r 3 ) := Xo 



'(ti,t 2 ,t 3 ) 



(6a) 




FIG. 1. Particle-hole scattering. 

(ii) In the pp-case we look at the scattering of two elec- 
trons with energies v' and lo — v'. Again the total energy 
of this process is u>. 



X^v(Tl,T 2 ,T 3 ) := X<TT'<T'a(ri,T2,T 3 ), 



(6b) 



which cover all six cases mentioned above. Eventually, 
using the crossing symmetry^, one can show that the 
quantity defined in Eq. (|6b[) can be obtained from the 
one given in Eq. (|6a[) by means of a mere frequency shift 
as it is explained in Appendix ID 21 For this reason we 
will commit ourselves to Eq. (pa|) and consider Eq. (|6bp 
only later when dealing explicitly with the spin-structure 
of the irreducible vertices. 

When switching to frequency space, it is convenient 
to define the Fourier transform of \ m two different 
ways, which we refer to as particle- hole (ph) and particle- 




(w - v) 



FIG. 2. Particle-particle scattering. 

Since in the full two-particle Green's function both pro- 
cesses are included, it is possible to express the \vv m 
terms of Xph and vice versa 



vv'uj _ uv'{u-u-v') 
Xpp^aa 1 Xph,acr' 



Xph 



vv'u vv' (u+v+v 1 ) 



(8) 



In the following, we will constrict ourselves to Xph = X 
and return to Xpp only when explicitly needed (all the 
definitions, results etc. of the following section apply 
also to Xpp)- 

In the case of an interacting system (U ^ 0), the sus- 
ceptibility x can be decomposed into two parts, in order 
to divide the bubble terms (independent propagation of 




A G„{ V ) 



GA"T 



FIG. 3. Diagrammatic representation of the generalized sus- 
ceptibility Xctct" i as defined in Eqs. (|7a|l and <(9j . In the inter- 
acting case Xctct'" i s naturally decomposed into a bubble term 
(Xo, see Eq. (fTO)) ') and vertex correction terms (F). 



the two particles) from the vertex corrections, as it is 
illustrated in Fig. [3] 



- ~f3G a {v)G a (v + u))5 UU '5 aa i 
G a {v)G„{p + u)FZ>G a ,{v')G ! ,,{v' + u). 



(9) 



The full vertex function F appearing on the r.h.s. of 
Eq. (jSJ) includes all possible vertex corrections, or in 
other words, all possible scattering events between the 
two propagating fermions, and can be hence interpreted 
in terms of the amplitude of a scattering process be- 
tween two quasi-particle a 33 ' 34 , at least in the Fermi-liquid 
regime, where the one-particle excitations are unambigu- 
ously defined. Eq. ([9]) can be also more compactly writ- 
ten in terms of the "one-particle" -like bubble part of x, 
defined as 



XT" 



-PG a (v)G a {v + uj)8 v 



(10) 



where the spin-indexes on the l.h.s. can be omitted by 
restricting oneself to the paramagnetic case 



VV UJ VV Ul 

Xaa' — AO "CTCT 



L \ * W!U T?VlV 2 U> V2V'0J 

aa ~fp 2^ X o o' Xo ■ I 11 ) 



Analogous definitions can be introduced for the particle- 
particle notation. 



B. Diagrammatics and mutual relations 

The full vertex-function F defined in Eq. ([9]) is the con- 
nected part of the complete four-point function. From a 
diagrammatic point of view F consists of all "fully con- 
nected" two-particle diagrams, i.e., all diagrams which 
are not separated into two parts. These diagrams, in 
turn, can be classified with respect to the way how they 
can be split into two parts by cutting two internal Green's 
function lines. 

(i) fully irreducible: Diagrams of F, which cannot be split 
into two parts by cutting two internal Green's function 
lines. They represent the two-particle "counter-part" of 
the self-energy diagrams at the one-particle level. 




FIG. 4. Schematic representation of a generic particle-hole 
reducible diagram contributing to the (full) scattering ampli- 
tude F. 



(ii) reducible: Diagrams of F, which can be split by cut- 
ting two fermionic lines. At the two-particle level, how- 
ever, the concept of reducibility is more articulated than 
at the one-particle level. In fact, there are more possibil- 
ities of cutting lines than in the one-particle case, and, 
therefore, the concept of reducibility has to be referred 
to a specific channel: This specifies in which way two 
of the four outer legs of a given diagram can be sepa- 
rated from the other two. Labeling the outer legs of the 
two-particle diagrams with 1, 2, 3, 4, it is clear that three 
different possibilities exist: If the outer legs 1 and 3 de- 
note outgoing particles (and 2 and 4 the incoming ones) 
than the diagrams where (13) can be separated from (24) 
are called particle-particle reducible, while the two other 
cases, i.e., (12) from (34) and (14) from (23), correspond 
to particle-hole longitudinal(p/i) and transverse (ph) re- 
ducible diagrams, respectively. One example for a (longi- 
tudinal) particle-hole reducible diagram is shown in Fig. 
|4] where (12) can be separated from (34) by cutting the 
internal lines a and b. 

It is worth recalling that each diagram is either fully 
irreducible or reducible in exactly one channel, i.e., there 
are no diagrams that are reducible in two or more 
channels^. As a consequence, the complete vertex func- 
tion F can be decomposed into four parts - a fully ir- 
reducible part A and the reducible contributions $ r in 
the three different channels (one particle-particle and two 
particle-hole) 



F = A + <I> T 



ph 



ph > 



(12) 



which have been written by now with a schematic nota- 
tion, omitting spin- and frequency arguments (they will 
be explicitly introduced in the next sub-section). 

Such a decomposition of F is known as parquet 



cquatio: 



,34,38 



and it is schematically illustrated in Fig. 



[5] with one low-order diagram shown for each of the four 
contributions. Note that the parquet equation represents 
just a "classification" of all connected two-particle dia- 
grams in four classes, and, therefore, does not imply in 
itself any kind of approximation. This is analog to the 
one-particle case, where all connected one-particle dia- 
grams can be divided into a set of reducible and irre- 
ducible ones (defining the self-energy). 

The full vertex function F appearing in the parquet 
equation can be calculated from the complete four-point 
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FIG. 5. Parquet equation. 



A 



all diagrams, e g : 





+ 



cduciblc in particle-particle 
channel 13-24 



ph 



+ 




reducible in particle- hole 
channel 12-34 




dble in particle-hole 
channel 14-23 



matrix element x ym Eq. © . In order to work with the 
parquet equation one needs additional relations connect- 
ing F and the reducible vertices $ r . This can be achieved 
by defining new quantities T r 



f = $ r + r r; 



= pp,ph,ph. 



(13) 



Since F contains all diagrams and $ r contains all the di- 
agrams which are reducible in the given channel r, T r 
is the set of all diagrams, which are irreducible in a 
given channel r. Since each diagram is either fully ir- 
reducible or reducible in a given channel, we have that 
r r = A + $ h +<$> jv jiJt^r. 

The T r vertices, in turn, can be calculated from F 
by means of an integral-equation, the so called Bethe- 
Salpcter equation: 



F = T r 



I 



T r GGF, 



(14) 



where the integral-symbol denotes an integra- 
tion/summation over all internal degrees of freedom 
(e.g.: frequencies, spin, ...). The interpretation of this 
equation is very simple: F is the sum of all connected 
diagrams which are irreducible in the given channel r 
(i.e., r r ) and the diagrams that are reducible in this 
channel (i.e., $ r ). The latter can be easily expressed 
by connecting the corresponding irreducible vertex T r 
to the full vertex function F, via two Green's function 
lines. Such a "decomposition" procedure, which avoids 
any possible double-counting of diagrams, is obviously 
not unique, as it can be performed independently for all 
channels. 

The two particle-hole channels are connected by the 
crossing relations^ 



i>(1234) = -1^(1432) 



(15) 



which corresponds to interchanging the two incoming 
particles. In contrast the particle-particle channel fulfills 
a crossing relation on its own, namely 



TABLE I. Approximations at different vertex levels. Vt r ) 
F, T r or A, respectively. 





static (V( r ) = U) 


local, dynamic 
(V(D=V ( ^ oc ) 


F 


2 nd -order perturbation 
theory 




Tr 


RPA, FLEX, 
pseudopotential parquet 


Moriyasque DrA 


A 


parquet approximation 


DrA 



r pp (1234) = -r pp (1432) 



(16) 



which is identical to the crossing relation for the full 
vertex F (cf. Appendix ID 2|) . 

In order to clarify the meaning of the different re- 
ducible and irreducible vertex functions (V( r ) = F, T r 
or A), defined in this section, it is important to discuss 
how some well-known approximation schemes correspond 
diagrammatically to different levels of approximation of 
the two-particle vertex functions. 

An overview over approximations adopted for the 
Hubbard model at different vertex levels is given in Tab. 
Ill 

The simplest schemes are obtained, obviously, by 
remaining at the "surface" of the two-particle diagram- 
matic complexity, i.e., when making approximations 
directly at the level of the full vertex function F. For 
instance, replacing F with the bare Hubbard interaction 
U leads directly to second order perturbation theory 
for the self-energy, as it can be easily seen by making 
such a replacement (F = U) in the Schwinger-Dyson 
equation of motion (cf. Eq. (f2"51 and Fig. IT61 below) . 
This method is used to approximate the self-energy of 
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the Hubbard model in the asymptotic cases of weak 
and strong coupling (see also, at half-filling, the iterated 
perturbation theory 8 -^) . 

Making a step further in the diagrammatics means to 
apply an approximation at the level of the irreducible 
vertices L r (third row of Tab. |T|. For example, one 
can calculate the full scattering amplitude F by simply 
replacing T p h (or T pp ) with the bare interaction U, 
which corresponds to the well-known random phase 
approximation (RPA)i£. Adopting this substitution 
for all three irreducible channels (and the fully irre- 
ducible one) leads to the Fluctuation Exchange (FLEX) 
approximatio n 41 ! 42 . If one chooses a different "effective" 
constant T r = C/Jf* f for each of the three channels, where 
the U^ff is determined by some additional condition, 
one ends up with the so-called pseudopotential parquet 
approximatio n 34 ' 43 . All these methods represent reason- 
able approximations in the case of small-intermediate 
[7, i.e., in the weak- intermediate coupling regime, and 
improve systematically the second order perturbation 
results. 

Remaining at the same level of the diagrammatics, 
a more complex approach, which aims to include non- 
perturbatively the physics of the MIT, is to replace T r 
by its purely local counterpart, instead of the lowest 
order (and frequency- independent) contribution U only. 
This is done in the ladder (or Moriyasque) version of 
the DLA2£. There, the reducible vertex $ r (typically 
in one or two of the different channels) is allowed 
to be non-local. Such an approximation is justified, 
when one of the channels dominates the physics (e.g., 
the spin channel for the half-filled Hubbard model 
at intermediate-to-strong coupling): In that case, by 
neglecting the interference between different channels, 
one constructs a non-local $ r from the local T r only in 
selected channels, while the remaining ladders can be 
considered as purely local quantities. 

Finally, going to the deepest level of the two-particle 
diagrammatics, one may apply approximations directly 
to the fully irreducible vertex A (fourth row in Tab. UJ). 
Again, for the Hubbard model the simplest approxima- 
tion of this class is obtained by replacing A with the bare 
Hubbard interaction U . This approach is called parquet 
approximatio n 34 ' 38 ! 44 . Similarly as for L r , however, one 
can also replace A with its purely local (but frequency de- 
pendent) counter-part, which corresponds to the DLA— . 
While the approximations at the level of A are usually 
very expensive computationally, and particular numeri- 
cal tricks 4 ^ have to be used, they may be necessary to 
capture the complicate physics of the Hubbard model 
in situations where none of the channels really domi- 
nates over the others (e.g., for the doped cases, which 
is relevant for the physics of the high-temperature su- 
perconducting cuprates). It is also worth recalling that 
the fully irreducible diagrams A have a very compact 



(y + uj) t (y : + u) t 




FIG. 6. Upper panel: lowest order (perturbative) diagrams 
for Fff, Lower panel: the same for Ffi 



structure, which cannot -per definition- include any lad- 
der diagrams. Therefore, for the case of the Hubbard 
model (where the naked interaction is completely local 
in space) it is reasonable to expect a weak spatial de- 
pendence, as it seems to be confirmed 4 ^ by a dynamical 
cluster approximatio n 46 ' 47 study for the two-dimensional 
Hubbard model. 

To conclude the diagrammatic classification of several 
known approximation schemes, it is worth to mention 
also the cases of DMFT— and DF— , which -in a strict 
sense- do not belong directly to any of the specific levels 
discussed above. In fact, DMFT is an exact theory in 
the limit d — > oo, and all local vertices (F, L r , and A) 
are included in its diagrammatics. However, the internal 
Green's function lines are also local in DMFT and non- 
local correlations are totally neglected in contrast with 
the methods discussed previously (from perturbation the- 
ory to DLA). On the other hand, DF does include spatial 
correlations beyond DMFT, via an expansion in a dual 
fermion space, defined via a Hubbard-Stratonovic decou- 
pling of the non-local degrees of freedom in Eq. ([1]). The 
coefficients of the DF expansion are given by the local 
generalized susceptibility of the AIM, which would cor- 
respond to a dynamical approximation at the level of the 
full vertex function F. However, a classification of DF in 
Table II would be not easy, as the local but dynamical 
vertex F represents the "naked interaction" in the dual 
fermion space. Hence, different diagrammatic degrees of 
accuracy are obtained by applying specific approxima- 
tions (perturbation theory/ladder/Parquet) for the dual 
fermions. While this makes a classification of the DF 
in Tab. Q] difficult at the moment, our discussion calls 
for future investigations of the correspondence between 
a given approximation in DF and the diagrammatics of 
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the real electrons. 



Spin-dependence: Definition of the different 
channels 



In this section, the spin-structure of the irreducible 
vertex in the three different channels is explicitly dis- 
cussed. As mentioned before, for the SU(2)-symmetric 
case there are three independent spin-combinations, i.e., 
ft, t-J- and tli see Eq^_©. On the level of the full ver- 
tex F the t^- and the fi-spin-combination are connected 
by the crossing relation Eq. (|D4bp given in Appendix 
ID 21 However, since at the level of T the crossing relation 
connects two different channels (sec Eq. (|15|l) at least in 
the particle- hole case we will for now consider the T r ^j 
as an independent quantity. Hence, we have three dif- 
ferent spin-combinations for each of the three channels, 
which would lead to nine different T's. Using crossing- 
and SU(2)-symmetry, however, one can show that only 
four of them are independent, which corresponds to the 
definition of the so called density (d), magnetic (m), sin- 
glet (s) and triplet (t) "channels" , given as follows 




FIG. 7. (color online). Vertex functions vs. the two fermionic 
frequencies v and v ': density part F% v w — U (left) and 
magnetic part F™' u + U (right) for U = 0.5 at half filling 
(/? = 26.0) at fixed cj; Upper panel: uj = 0, Lower panel: 
w = 40f. 
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where for the two particle-particle channels the 
corresponding particle-particle frequency definition is 
adopted. The first two quantities correspond to the lon- 
gitudinal particle-hole channel, i.e., T p h while the third 
and the fourth definition are linear combinations of ver- 
tex functions irreducible in the particle-particle channel. 
It is not necessary to explicitly consider the transverse 
particle-hole channel T^, since it is connected with the 
longitudinal one by means of the crossing symmetry Eq. 

G9. 

The same definitions are valid at the level of F and A 
as well. However, since neither F nor A can be divided 
into different channels, only two of them are actually in- 
dependent. 

A more detailed discussion of the different irreducible 
channels can be found in Appendix |B] and in Rcf. l34l . 



III. DMFT RESULTS 



Noteworthily, our DMFT(ED) results for reducible and 
irreducible local vertices have been successfully compared 
with corresponding results obtained with a Hirsch-Fye 
quantum Monte Carlo algorithm^ as impurity solver, in 
a slightly higher temperature regime than that consid- 
ered here. 

Following the thread underlying the discussion of Tab. 

II at the end of Sec. II B, we will start analyzing in 
the next subsection the most conventional (and easiest 
to compute) among the vertex functions, i.e., the full 
vertex F. Subsequently, in Sec. Ill B we will make a 
step deeper in the diagrammatics, presenting our DMFT 
results for the irreducible vertices in one specific channel 
T r , and, finally, in Sec. Ill C, results for the most fun- 
damental block of the two-particle diagrammatics, the 
fully irreducible vertex function A, will be presented and 
discussed. 

In all cases, the frequency structure of the local vertices 
will be first examined at small values of U (e.g., U = 0.5), 
which allows for a direct comparison with perturbation 
theory. Deviation from the perturbation theory predic- 
tions will be also discussed, especially at the end of Sec. 

III C, where their effects on more conventional physical 
and thermodynamical quantities will be eventually ad- 
dressed. 



We present in this section our new DMFT results for 
all local two-particle vertex- functions, i.e., F (full ver- 
tex), T r (irreducible in channel r) and A (fully irreducible 
vertex) of the half-filled Hubbard model on a cubic lat- 
tice. The frequency-dependent local vertex functions 
have been obtained by solving the AIM associated to the 
DMFT solution by means of exact diagonalization (ED). 



A. Full vertex functions 

The full vertex F contains all connected diagrams with 
two particles coming in and two particles going out. 
In Fig. [S] the lowest order diagrams for the two pos- 
sible spin combinations are shown in the particle-hole 
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frequency convention (the corresponding results in the 
particle-particle notation can be simply obtained via the 
transformation w — > lj—v—u'). Wc recall, moreover, that 
on the level of F, the singlet- and the triplet-channel are 
just linear combinations of Fd and F m . 

In terms of Green's functions the lowest order contri- 
butions for F (but it is the very same also for T r and A) 
read as follows 



Vi 
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for the fi-case. The lowest order contribution for the 
four different channels defined in Eq. (fT7j) -(|20 |) hence are 
given by 
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In the upper panel of Fig. [7] the full vertex functions 
F in the density (Fd — F^ + F^±) and magnetic 
(F m = Fff — Ffi) channel calculated by means of 
DMFT are plotted for the case u> = 0. The a;-axis 
corresponds to v while the y-axis is assigned to v' . Note 
that, for the sake of readability of the figure, instead 
of the absolute values of the Matsubara-frequency 
just the corresponding indexes are given. The vertex 
functions F are calculated for U = 0.5 at half filling, 
at a temperature value (j3 = 26.0) close to the critical 
end-point of the metal-insulator transition in DMFT—. 
It should be recalled that for the half-filled system 
all vertex functions are purely real (see Eq. (|D19|) in 
Appendix ID 4j) . Furthermore, here as in the following, 
the (constant) contribution of the first order diagram, 
namely the Hubbard U, is subtracted in order to better 
highlight the frequency structure of the two-particle 
vertices, beyond the standard lowest order perturbativc 
results. 



One can now trace the different features of the 
two-dimensional plot of F back to different types of 
diagrams. 

First of all, let us note that a constant background is 
still present, despite the subtraction of the lowest order 
term. This constant background stems from higher order 
diagrams that are independent of v and v 1 . An example 
in second order perturbation theory is given in Fig. [6] 
The left diagram in the upper panel ('Pi) lias no v- or 
^'-dependence, as it also follows from Eq. (|21ap . The 
same holds also for diagrams of higher order with all 
possible vertex corrections inside the bubble of V\ : The 
sum of all diagrams of this family yields the constant 
background observed in the upper panel of Fig. [7] 
Secondly, the evident structure along the main diagonal 
(i.e., the region around the line v = v') stems from 
diagrams like the second ones (V2, "Pa) in the upper or 
lower panel of Fig. [S] (sec also Eqs. (|21b|) and (|22b|> ). 
which describe (at the order considered) scattering 
processes reducible in the particle-hole channel. More 
specifically, these diagrams, as well as similar diagrams 
of the same type but with vertex corrections included, 
depend only on (y — v'), which means that they give 
a constant contribution along the lines v — v 1 = const. 
The largest contribution, however, is expected for the 
case v — v' = when the scattering between the particle 
and the hole occurs at the Fermi surface. One can easily 
identify these structures in Fig. [7] For the density-case 
one obviously has to add the diagrams of the ft - an d 
the fi-channel which leads to twice the contribution 
of such diagrams in second order perturbation theory. 
For the magnetic vertex, instead, these second-order 
contributions cancel exactly each other, and only higher 
order contributions to this diagonal line remain, which 
explains the difference between the two channels. 

Furthermore, one also observes an enhanced scatter- 
ing rate along the secondary diagonal v' = — v (for the 
magnetic channel its value is reduced due to cancellation 
with the constant background) . The origin of this struc- 
ture stems from diagrams like V5 in the lower panel of 
Fig. [6] (see also Eq. ()22c|) ). which build up scattering 
processes reducible in the particle-particle channel. In 
fact, such diagrams (with and without vertex corrections 
in the bubble) describe the scattering of two particles 
with energies (y + w) and v 1 . Hence, the corresponding 
scattering amplitude is enhanced for total energies at the 
Fermi level, i.e., for v' = — v — oj. If oj = this yields the 
secondary diagonal in the plots shown in the upper panel 
of Fig. [71 However, for a finite uj this line is expected to 
be shifted to v' = —v — u. This behavior is shown for 
case of the twentieth bosonic Matsubara frequency, i.e., 
for w = J(2 x 20) = 40|, in the lower panel of Fig. [7j 
The main diagonal remains unchanged, as it stems from 
w-independent diagrams, while the secondary diagonal is 
shifted compared to the upper panel. 
Finally, one can also note a "cross-structure" in the upper 



9 




FIG. 8. (color online). Irreducible vertices in different channels. Upper panel: Y v d " w -U,F% u + U, u - 2U, V v t v u for 
U — 0.5, /3 = 26.0 and uj = at half-filling. For singlet- and triplet-channel particle-particle notation was adopted. Lower 
panel: one-dimensional snapshot of the same vertex functions for v = 5 (fixed) and two different values of cj (ui = and 



20f) 



(v 2 + v' - u t ) t 





FIG. 10. (color online). Irreducible particle-particle vertices 
in particle-hole notation: < ~ v+u + "'— 2(7 (left panel) and 
T"" '" + " + "' (right panel) for the same parameters as in Fig. 

m 



FIG. 9. Third order (perturbative) diagram for F-pf 



panel of Fig. [7J (which exists, but is much weaker, in the 
density channel too), i.e., one observes an enhanced scat- 
tering amplitude compared to the constant background 
along the lines v — and v 1 — 0. In order to explain 
the origin of these structures one has to go at least to 
third order perturbation theory: The contribution of the 
diagram shown in Fig. [9] reads as 

Pe = ^ J2 G{ Vl )G( Vl +u:)G(v 2 )G( V2 + V ' (27) 

One sees that it is independent of v and therefore it 
gives a constant contribution along the horizontal line 
v' = const, in Fig. [71 with a maximum for v' ~ 0, as 
only in this situation one has the possibility to have all 
Green's functions appearing in Eq. (|27[) simultaneously 
at the Fermi level. In complete analogy one can con- 
struct diagrams that do not depend on v' (in this case 



the vertical bubble should be on the left side). These 
result in a maximum at u = 0, which completes the 
explanation of the "cross-structure" observed in F. 



B. Irreducible vertices in one selected channel 

At the level of irreducible vertices one has necessarily 
to consider four independent quantities: the density and 
the magnetic vertex correspond to the two possible spin- 
combinations in the longitudinal channel (T p h) while the 
singlet and the triplet vertex are linear combinations of 
the two different spin-directions in the particle-particle 
channel (r pp ). The transverse channel is not indepen- 
dent since it can be obtained from the longitudinal one 
by means of the crossing symmetry (see previous section 
and Appendix). 

We start with the discussion of the two-dimensional 
density-plots for the four channels for U = 0.5 and 
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io = (Fig. [8]). It is important to recall that for the 
two particle-particle channels, i.e., the singlet- and the 
triplet-channel, the particle-particle notation is adopted. 
For the density- and the magnetic channel (see upper 
panel in Fig. [8]) one identifies the main and the secondary 
diagonal as it was the case for the full vertex function. 
However, the constant background and the "cross" in the 
center are missing. In fact, such features originate from 
diagrams like the first one in the upper panel of Fig. [6] 
which are reducible in the longitudinal channel and there- 
fore do not contribute to and r m . 
For the particle-particle channel one would expect again 
a constant background (as well as the "cross" structure) 
but no secondary diagonal since the last diagram in the 
lower panel of Fig. [6] does not contribute (it is particle- 
particle reducible). This is, however, not true in the 
particle-particle notation, i.e., oj — >• u — v — v' ': The first 
diagram V\ in the upper panel of Fig. in fact, depends 
on uj — v — v' (instead of being independent of v and v' 
at all) and therefore yields a constant contribution along 
the lines ui = v + v' . For the case ui = this contribu- 
tion reaches a maximum yielding the secondary diagonal 
structure, as it appears in the density-plots for Y s and 

r t . 

On the other hand, if one uses the particle-hole nota- 
tion also for the singlet- and the triplet vertex one gets 
the constant background as one can see in Fig. [TO] Let 
us mention another interesting feature characterizing the 
triplet vertex T t in the particle-hole notation: The triplet 
vertex T t coincides with the '("'f-vertex. Hence, in the 
particle-hole notation, it describes the effective interac- 
tion between two electrons with spin f and energies v + ui 
and v' ', respectively (these are the energies associated 
with the two annihilation operators in the particle-hole 
notation). However, for v' = v + u> both electrons would 
be in the same state, which is forbidden by the Pauli- 
principlc. Therefore the triplet- vertex is expected to be 
strongly suppressed along this line, as it can be actually 
observed for ui = in Fig. [TUJ (right panel). 
In the lower panel of Fig. [8] one-dimensional slices of 
the four irreducible vertex functions are shown: v is kept 
fixed to the first fermionic Matsubara frequency in that 
case, and T is plotted for two different values of a; as a 
function of v' . One observes a good agreement with de- 
viations of order U 3 ~ 0.1. This is to be expected since 
third order diagrams have not been considered in the 
perturbation expansion. At U = 1.0 quantitative devia- 
tions from perturbation theory results become gradually 
visible in the low-frequency (v, v') region, see Fig. Ill| 
especially for the density channel. However, the main 
structures of the Y vertices radically changes w.r.t. per- 
turbation theory only at U ~ 1.5. 



C. Fully irreducible vertices 

In this subsection we present results for the fully ir- 
reducible vertex A. The formulas used for the actual 



calculations are given in Appendix [C] As mentioned be- 
fore, the fully irreducible vertex is the most fundamental 
"brick" among the two-particle vertex functions, repre- 
senting the diagrammatic analog of the self-energy, but at 
a two-particle level. Hence, approximations based on this 
level of the diagrammatics, such as the parquet approxi- 
mation or the DrA are extremely appealing from a theo- 
retical point of view. At the same time, the calculation of 
fully irreducible vertex functions is quite challenging, so 
that the few calculation o 38 i 44 ' 45 based on approximations 
for A, simply replace the latter with its lowest-order con- 
tribution (U). Motivated by the lack of studies on the 
frequency dependence of the fully irreducible local ver- 
tices, even at the level of perturbation theory, we will 
present our numerical and analytical results with more 
details than in the previous subsections and we also ex- 
plicitly consider the effects of the frequency dependence 
of A in selected physical and thermodynamical quantities 
as a function of the Hubbard interaction U. 

By definition no channel-dependence of the fully irre- 
ducible vertex function A can exist, since it is irreducible 
in all channels. Hence, as in Sec. Ill A, here we also 
restrict ourselves to the DMFT result for the density 
and the magnetic vertices, which represent the two pos- 
sible spin combinations. Diagrammatically, the lowest 
order contribution to the fully irreducible vertex is the 
bare Hubbard interaction U (Fig. [6]). The next terms 
in the perturbation expansion are already of 4 th order: 
These diagrams have the form of an envelope, and, hence, 
arc usually referred to as "envelope" -diagrams. The 
envelope-diagrams for the f'f- and the f^-case are shown 
in Figs. (fl~2j) and (fT3|) . respectively. 

Let us just mention one interesting feature for the ft 
diagrams, which is relevant for the particle-hole symmet- 
ric case discussed here: At half filling the first and the 
second (as well as the third and the fourth) diagram be- 
come exactly the same, i.e., one can take only the first 
and the third diagram and assign a factor 2 to them. 
This happens because all these diagrams differ only for 
the direction of the propagators in their closed fermion 
loops containing four internal electron ^-lines. This is 
analog to Furry's theorem^ in quantum electrodynam- 
ics, which states that -as a results of the electron-positron 
symmetry- a closed fermionic loop containing an odd 
number of Fermions always vanishes, while a closed loop 
with an even number of fermions gets a factor 2. 
In Fig. [TU eventually, our DMFT results for k u / u and 
A^ w are compared with the the ^-contributions from 
perturbation theory, given by the envelope diagrams in 
Figs. [T2l and [T51 respectively. Algebraically, the contri- 
bution stemming from such a diagram is 

A e ™ = (±)^ Y, G{ Vl )G{v 2 )G{v 3 )G{ Vi )G{ V5 )G{ V6 ), 

(28) 

where 1/4, v 5 and v e are functions of f\, v<x and 1/3 
(reported in Figs. [T^] and [T5|) rather than independent 
summation variables. We recall, that the lowest order 
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diagram, which is simply given by the bare Hubbard 
interaction U (first diagram in the second panel of Fig. 
[5]), is subtracted in both cases, i.e., only deviations from 
this constant contribution are plotted. From Fig. [14] one 
can see that the structure of the exact results for A^ m " 
resembles very much that of the envelope diagram. This 
is expected for a relatively small U = 0.5 and it is also 
demonstrated in Fig. 1151 There a one-dimensional slice 
of A is plotted (lu and v are fixed) in comparison with 
perturbation theory, i.e., the envelope-diagrams: The 
deviations from the constant term ±U for k" d v u and 
A™ u , respectively, are of the order U 4 ~ 10~ 2 — 1CP 3 
which is perfectly consistent with our numerical data. 
Let us furthermore mention that our numerical results 
demonstrate that the high-frequency asymptotics of the 
fully irreducible vertex A reduces to the lowest order 
perturbative contribution (U), which is consistent with 
the analytical derivation of the high-frequency behavior 
of T in Ref. [H. 

At the end of this section, we establish a connection 
between the results for two-particle quantities, we have 
presented so far, and the more familiar results at the one- 
particle level, i.e., those for the self-energy of the system. 
In a second step, the connection with some selected physi- 
cal quantities, which are typically analyzed in the context 
of the Hubbard model, will also be illustrated. 

As for the self-energy, this goal can be easily achieved 
by exploiting the Heisenberg (or Schwinger-Dyson) equa- 
tion of motion 

EM = ^-^X^W)G(z/+ W )G(«/+ W ). (29) 

whose diagrammatic representation is given in Fig. 1161 

Inserting the parquet Eq. (fT2"j) in Eq. (|2"9j) , i.e., 
splitting up u into a fully irreducible and the 
three reducible parts allows us to identify four different 
contributions to the self-energy stemming from the 
irreducible and the reducible part of the full vertex 
function F. Specifically, in Fig. (fT7|) we compare the 
DMFT self-energy E(f) with its contribution which 
originate from the fully irreducible vertex A only and 
from its lowest order contribution (U), respectively. For 
this purpose we used Eq. (|29|) and replaced F^ " by 
A'fY w and U, respectively (where the second case simply 



yields diagrammatic contributions similar to those of the 
2 nd -order perturbation theory). 

For the relatively small U = 0.5 (left panel) there 
is no visible difference between the calculations of the 
self-energy with the full A and U. This is to be expected 
here, since in the perturbative regime the relative 
difference between the full A and U is extremely small, 
as we have already noticed in Fig. [151 

On the other hand, for a larger value of the Hubbard 
interaction (U = 1.5, upper right panel in Fig. 1 1T[) one 
can see that, while the major part of the self-energy is still 
coming from the fully irreducible part of F, the frequency 
dependence of A becomes essential for the calculation of 
the self-energy. Hence, setting A = U, as it is done, e.g., 
in the parquet approximation, would yield results quite 
far from the correct structure of the one-particle local 
self-energy, which appears to be determined to a major 
extent by frequency-dependent high-order terms of A. 

Very similar conclusions can be drawn by analyzing the 
contribution of A to the value of the double occupancy 
rifni = as a function of U, which is 

shown in the lower-left panel of Fig. 1171 By compar- 
ing the results obtained with the full DMFT self-energy, 
with those obtained considering A only, we observe that 
also in this case the irreducible vertex gives a signif- 
icant contribution to the well-known reduction of the 
double-occupancy value w.r.t. its not interacting value 
of n-j-n^ = n-|- x nj, = 0.25 with increasing U. Also in this 
case, however, for U > 1 the results calculated with the 
approximation A = U deteriorate very quickly, so that at 
U ~ 1.4 a very incorrect estimate for n^n^ would be ob- 
tained by neglecting the frequency dependence of A. The 
situation appears more articulated, however, when ana- 
lyzing the case of two-particle local response functions, 
such as the density Xd" 1 an d magnetic Xm C local suscep- 
tibilities at zero bosonic frequency (uj = 0). Such ther- 
modynamic quantities contain a very important piece of 
information for the physics of the Hubbard model: Ap- 
proaching the MIT is marked by a constant enhancement 
of Xm C w ith increasing U. In fact, a divergence of Xm C 
actually signalizes the transition line, as it corresponds 
to the formation of a stable local magnetic moment in 
the Mott phase. At the same time, the reduced mo- 
bility of the electrons with increasing value of the local 
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FIG. 12. (color online). [^-contributions to the perturbative 
expansion of the fully irreducible vertex A"" " ( "envelope" - 
diagrams) in particle-hole notation. 



Coulomb interaction U is mirrored by a gradual suppres- 
sion of the local charge fluctuations, and, hence, by a 
monotonous decrease of Xd° c > with U. Such trends are 
naturally well captured by our DMFT calculations, per- 
formed via a summation of both Matsubara fermionic 
frequencies v, v' of the generalized susceptibility \ vv 
defined as in Eq. (|9|). Following the same procedure de- 
scribed above, we have singled out the contribution to 
\ l £ c and x l m originated by the fully irreducible vertex A 
its lowest-order term (U). While only limited informa- 
tion can be extracted from the Xd° C ' as ^ is becoming very 
small in the non-perturbative region, by analyzing the 




V± = V + V + UJ — V\ + V2 — U3 

Vi t u 5 = V + UJ + v 2 — v 3 



FIG. 13. (color online). L^-contributions to the perturbative 
expansion for the fully irreducible vertex AST " ( "envelope" - 
diagrams) in particle-hole notation. 
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FIG. 14. (color online). Upper panel: density plots of A™ " — 
U (left) and A%'" + U (right) for U = 0.5, /3 = 26.0 and uj = 
at half filling; Lower panel: 4 th -order perturbation theory 
results ("envelope" diagrams in Figs. H2land [T3l respectively). 



data for \m some relevant difference with the previous 
cases can be noted. The contribution to Xm C stemming 
from the irreducible vertices A and reducible diagrams 
in F are comparable. The latter contributions appear to 
become the predominant ones in the region U > 1 where 
a stronger enhancement of \ l m is observed. We also note 
here that the relative error made by replacing A = U is 
naturally increasing with U but remains weaker than in 
the previous cases. 
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FIG. 15. (color online). K" d " u -U (left) and A™ " + U (right) 
for U — 0.5 (f3 = 26.0) for selected one-dimensional snapshot 
at fixed oj = or ui — 205, J-* = 5, as a function of v . 
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FIG. 16. Schwinger-Dyson equation of motion. 



DMFT RESULTS FOR THE ATTRACTIVE 
MODEL 



Exploiting the particle-hole symmetry of the half-filled 
Hubbard model, it is possible to perform an exact map- 
ping between the repulsive (U > 0) and the attractive 
(U < 0) Hubbard model. On the local level, which we 
are interested in, this mapping is performed by a unitary 
operator—, which transforms particles with a given spin 
(e.g., the J,-spins) into holes and vice versa, i.e., an oper- 
ator cj becomes (and vice versa) under this transfor- 
mation. Particles (holes) with the other spin-direction 
(i.e., f m this case) are invariant. The Hamiltonian is 
also unchanged except for a change in sign of the Hub- 
bard interaction parameter U . The same holds also for 
the purely local model, i.e., the AIM associated with the 
DMFT solution (see also Appendix IDT)) . 
One of the consequences of this symmetry is e.g., that the 
Green's function and the self-energy for systems which 
differ only in the sign of the Hubbard interaction U are 
identical. Furthermore, for the local model, one can show 
that these one-particle functions are purely imaginary. 
At the two-particle level the situation is logically more 
complicated. A detailed calculation (sec ID 4[) shows that 
the ^^-susceptibility calculated with U < is mapped 
onto to the magnetic susceptibility for U > (Eq. 
(|D23|) ). Physically, this can be understood as follows: 
Fluctuations of the x— and y— spin component (i.e., the 
order parameter of the antiferromagnetic phase transi- 
tion) at positive U are equivalent to fluctuations of the 
"cooper-pair density" c|c^ (i.e., to the superconducting 
order parameter) at negative U . For the lattice model 
this means that for an (antiferro)magnctic instability at 
a given point (U > 0,T) in the phase diagram there ex- 
ists a superconducting instability in the corresponding 



FIG. 17. (color online). Upper panel: DMFT self-energy (red 
circles) compared to the contributions stemming from A (blue 
squares) and from U (light-blue triangles) only, respectively, 
for U — 0.5 (left) and U = 1.5 (right); Lower panel: double 
occupancy (left) and susceptibilities (right). The error bars 
refer to the finite frequency range adopted for the fermionic 
frequency summations over v, v'—. 
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FIG. 18. (color online). Upper diagram: reducible in particle- 
particle channel, Lower diagram: reducible in the transverse 
particle-hole channel, — U denotes the attractive and +U the 
repulsive model 



attractive model at (—U,T). 

While a complete algebraic derivation is given in Ap- 
pendix [D4l we provide here a diagrammatic illustration 
of the relationship between the two channels (i.e., the 
magnetic and the particle-particle fl) ■ We start with an 
arbitrary fi diagram which is reducible in the particle- 
particle channel, i.e., it contributes to ^pp^(-u^ (see 
upper panel of Fig. [T5|) . The ^-Green's functions (plot- 
ted in red) are reversed under the particle hole transfer- 
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FIG. 20. (color online). T™' u for U = +0.5 (left) and 

V PP,n = K r r'" +rr' w ) for U = -0.5 (right), both for 
P = 26.0 and w = 



which state that the spin-channel is mapped onto the 
particle-particle t|-channcl (plus an additional frequency 
shift), we compare these two channels in Fig. [20] Per- 
forming the additional transformation v' — > (u> = 0) — v 1 
in the plot for r^J^ (i.e., reflecting the plot along the 
x = !/-axis) one sees that the two plots are identical, 
which illustrates numerically very well the attractive- 
repulsive mapping of the half-filled Hubbard model at 
the two-particle level. 



V. CONCLUSIONS 



mation. 

Naturally the corresponding frequency arguments also 
change their sign. The diagram we obtain after the 
particle-hole transformation is a diagram which is re- 
ducible in the transverse (=vertical) particle-hole chan- 
nel. Since this relation holds for all reducible diagrams 
we can formulate the following equation for the $'s 



(30) 



(The minus-sign stems from the exchange of the two 
fermions.) Furthermore SU(2)-symmetry states that 



lation in Eq. (|30|) yields 



(see Eq. (|D6|) V Using this re- 



v{y) — v' ) (— uj) 



(31) 



Finally the additional transformation v' — > uj — v' (sec 
Appendix ID 41 and [B| gives the mapping. Evidently this 
relation holds for the T's as well. This result is in com- 
plete agreement with Eq. (|D23[) since the T's are calcu- 
lated from the x's by an inversion. 

Results for the T's in the four different channels are 
shown in Fig. [19] for the case U = —0.5. Comparing 
it with Fig. [8] i.e., the T's for the corresponding repul- 
sive case U = +0.5, one observes that the triplet-channel 
" is unchanged. This is expected because the triplet 
channel is identical to the particle-particle channel 
^pp ft' an< ^ ^ e t _crea tion- and annihilation-operators are 
not affected by the particle-hole transformation. 
Furthermore, following Eqs. (|D23]l and (|D24| . 



The problem of electronic correlation still represents a 
big challenge in the field of condensed matter physics, and 
the theoretical analysis is very hard even for the easiest 
model of electronic correlation, i.e., the Hubbard model. 
In this respect, DMFT has represented an huge step for- 
ward, as it is able, e.g., to describe non-perturbativcly the 
opening of the Mott-Hubbard gap in the phase-diagram 
of the half-filled Hubbard model. The study of the quan- 
tum local fluctuations in DMFT, however, has been often 
performed at the one-particle level (self-energy, spectral 
function), or for quantities that can be directly extracted 
from the self-energy calculations (e.g., the optical con- 
ductivity). A precise knowledge and understanding of 
the two-particle local vertex functions is, however, crucial 
both for computing momentum-dependent two-particle 
spectral properties at the DMFT level and for diagram- 
matic extensions of DMFT, aimed at the inclusion of 
long-range spatial correlations, such as the DF and DrA. 

In this work, we have analyzed systematically the prop- 
erties and the frequency structures of the two-particle lo- 
cal vertex functions by means of DMFT, applied to the 
half-filled Hubbard model on a cubic lattice. The analy- 
sis of the two-particle vertices was not limited to the full 
scattering rate amplitude F or to the irreducible vertices 
in a given channel T r , which have been compared and 
interpreted in terms of the- corresponding perturbation 
theory results and of the mapping onto the attractive 
Hubbard model. On the contrary, for the first time -to 
the best of our knowledge- we have also presented DMFT 
results for the fully irreducible local vertices A. Specif- 
ically, it has been shown how the frequency dependent 
part of the fully irreducible vertex function, crucially af- 



fects the one-particle self-energy at intermediate values of 
the Hubbard interaction. This would imply that approx- 
imations at the level of the fully irreducible vertex, which 
rely on the parquet formalism, should necessarily include 
its frequency dependence, as, e.g., in the DrA, in order 
to go beyond the weakly-correlated regime. Whether, 
and to what extent, it is possible to neglect the spatial 
(or momentum) dependence of the fully irreducible ver- 
tices of the Hubbard model remains to be investigated, 
though a specific set of data obtained in dynamical clus- 
ter approximation for a two-dimensional Hubbard model 
appears rather promising in this respect. 
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FIG. 21. Domain of definition for G(ti,T2) 



Appendix A: (Imaginary) Time translational 
invariance - Boundary conditions 

We summarize here the so called Kubo-Martin- 
Schwinger (KMS) boundary conditions for the n- 
particle Green's function, which follow from the time- 
translational invariance and from the cyclic property of 
the trace (see Ref. [H[ and Ref. [3l|). 

We start from the n-particle Green's function defined 
in Eq. (J3j> omitting the spin-indexes for this section, 
since the subsequent considerations are independent of 
the spin. Note also that the results discussed here are 
valid for models with arbitrary degrees of freedom (like 
spin, lattice site, etc.) and not only for the local AIM 
Hamiltonian: The only requirement is time-translational 
invariance of the system, i.e., that the Hamiltonian T-L is 
independent of r. 

Assuming that t\ is the largest and T2„ is the small- 
est time argument of the n-particle Green's function 
G„(ti, . . . , Tin) one gets the following conditional for the 
In time- variables: 

T 2n + P > n > . . . > n > . . . > t 2 „, (Al) 

i.e., all time- arguments have to be located within an in- 
terval of the length /?. Otherwise the term e~^ +T2n ~ Tx ' H 
in the definition of G n would lead to exponentially in- 
creasing contributions with growing eigenvalues E n of 
the system, and the trace occurring in Eq. ([3]) diverges. 
On the contrary, if condition (|A1|I is fulfilled, the above- 
mentioned exponential factor suppresses the contribution 
to the trace for large eigenvalues E n of H and hence, 
the trace converges and the n-particle Green's function is 
well defined. As an example, the domain of definition for 
the one-particle Green's function Gi_(t\,T2) = G(ti,T2) 
is shown in Fig. (region between the two solid diago- 
nal lines). 

Due to the time-invariance of the Hamiltonian the n- 
particle Green's function G n does not depend on all 2n 
times explicitly but rather on time-differences, e.g., of 



the form r, — r 2n , yielding 

G„(ri, . . . , T 2n ) = G„(ti — T2n, ■ ■ • , T2n-1 — T2„, 0). (A2) 

As a result, the one-particle Green's function is constant 
along diagonals of the form r 2 = T\ + a, a G [— /3, /3] in 

Fig. HQ 

Furthermore, the cyclic property of the trace leads to 
anti-periodicity of the n-particle Green's function which 
reads 

G„(ti, . . . ,T 2 „) = -G„(ri -/?,.. • ,T 2n ), (A3) 

if we assume that 7*1 > ... > r 2n > t\ ~ (3. 
All imaginary time- variables can be restricted to the in- 
terval [0, (3], since the value of G n for all other combina- 
tions of time-arguments (that are allowed according to 
Eq. IA1[) can be constructed by means of Eqs. (|A2|I and 

Considering the anti-periodicity condition (|A3|) one can 
express the n-particle Green's function G n as a Fourier- 
expansion 

Gn(ri, . ■ • ,r 2 „) = -1 e- i{viTl+ ~- v * nT * n) G n {v u u 2n ), 

G(p 1 ,...,u 2n )= f dr x ... f ' dr 2n e^+--^ T ^x 
Jo Jo 

x G n (n, . . ., r 2n ), 

(A4) 

where Vi = ^(2n^ + 1) are fermionic Matsubara- 
frequencics. The calculation of the Fourier-coefficients 
can be simplified by means of the following considera- 
tions. One uses Eq. (| A2|) and performs the substitutions 
t, = t[ — T2„, i = 1, . . . , 2n— 1. Next, one can shift the in- 
tegration intervals of r{, . . . , T 2n _ 1 from [— r 2 „, ft — r 2n ] to 
[0,/3] due to the anti-periodicity condition (|A3[) . Hence, 
the r 2n integration in Eq. (|A4|) can be performed analyti- 
cally and leads (beside a factor f3) to energy conservation 
v\ — v 2 + . . . + V2 n -\ — v 2n = 0. Therefore, it is sufficient 
to consider a (2n — l)-frequency object 

fP fl 3 

Gn(vi, . . . ,V2n-l) = \ &T\ . . . / d,T 2n -l X 

Jo Jo 
xe i ^+--»*«-^-^G n (T 1 , . . . ,t 2 „_i, 0), 

(A5) 

related to the full 2n-frcqucncy Green's function via 

G n {v\, . . . , V2n) = /3S(v 1 +...+v 2 „_ 1 )v 2n G n (vi, . . . , V 2n -l)- 

(A6) 
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FIG. 22. Bethe-Salpeter equations in the longitudinal chan- 
nel. 



Appendix B: Spin-diagonalization 

In this appendix, we summarize the spin-dependence of 
the three (ir)reducible channels (i.e., ph,ph and pp) and 
give a derivation of the corresponding Bcthc-Salpctcr- 
equations for the SU(2)-symmetric case. The formalism 
presented here is similar to that of Ref . 34 1 . 



a. The longitudinal (horizontal) channel Y p h 

We start with the Bethe-Salpeter equations for the 
three independent spin-combinations T p h t ff, Yph^i and 
T h tt- Diagrammatically they take the form shown in 
Fig. [2H Algebraically they read as 

Fn'- = r^ t +~ E ^M^M^+^Kt;'" (Bla) 

Frr = ntn+- p E k^g^g^+^f^ (Bib) 

It is easy to verify the plus-sign in front of the sec- 
ond summand on the right hand side of these equations 
by comparison with 2 nd -order perturbation theory: The 
corresponding perturbative contribution shown in Fig. 
El upper left diagram, exhibits a plus-sign (see also Eq. 

(25)). 

One can see that Eqs. (|BTa[) and (|BTb[) are coupled, 
while Eq. (|Blcj) contains only Anyway, we will 

postpone the calculation of this vertex function to the 
transversal particle-hole case since r , jt is related to 
^pK fj. by the crossing relation Eq. (fT5|) which reads as 



pl^f UJ 

ph.Ji 



-r 



v{y-\-uj){i>' — v) 



(B2) 



for this specific case. 

By hands of SU(2)-symmetry the two other equations 
can be decoupled analytically considering the sum and 
the difference of Eqs. (|Blap and (|Blb[) . respectively: 



T?VV UJ 

d(ensity) 



t?w uj 



T?W UJ 



r m(agnetic) * — ^tt ti 1 



(B3a) 
(B3b) 



which correspond to Eqs. <jTTJ) for the T's. The two 
decoupled equations for the density and magnetic channel 
are 

F v / U = T? u + \ Y, V 7 l "G( Vl )G{ Vl + u)F v /»'", 

(B4a) 

(B4b) 

These equations can be solved for the T's by an inversion 
of the matrix (1 + ±GGFy u ' u in the zV-space, i.e., 



IT™ = E C: K 1 + -pGGF d , m )^} ™ 



(B5) 



Considering the definition of x m Eq. ([9]) one can write 
the quantity which is inverted as xYm I Xq v " ■ 
For the sake of completeness, we want to rewrite this 
equation into the form which was used for extracting the 
T's shown in this paper. Defining xY u an d Xm " an d 
combining Eq. ([9]) with Eqs. (|B4|) one finds the corre- 
sponding Bcthe-Salpeter-equations for the x's: 



VV UJ _ VV UJ 

Xd.m AO 



02 2—4 *° d,m Xd,m ' V-° U J 

V1V2 



Solving these equations for T" d u w and " yields 

fiXdl-Xo'r'"- (B7) 



-pVV UJ 



b. The transverse (vertical) channel 



The Bethe-Salpeter equations for the three different 
spin-combinations shown diagrammatically in Fig. [53] 
read as 



x F- 



KwiK^+uiKu-ui) 



o-lt 



Sa) 
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FIG. 23. Bethe-Salpeter equations in the transverse channel 
with v = v + u\, v = v + uji , Q = uj — wi. 



(B8b) 



x R 



WlCTl 

ji(f+«l)(»'+ui)(u-ui) 



x F± 



(B8c) 

As in the longitudinal channel the minus-sign in front of 
the reducible part of these equations can be inferred from 
comparison with 2 nd order perturbation-theory (vertical 



diagrams in Fig. [6] as well as Eqs. ()21b[) and (|22a|l ). 
One can see that Eqs. (|B8aj) and (|B8c[) arc not inde- 
pendent, i.e., in the transverse channel the ft - an d the 
ti- vertex are coupled in the same way as it was the case 
for T p h^ and Tph,fl m the longitudinal channel (see Eqs. 
(|Bla|) and (|Blb|) ) . This is not surprising since these func- 
tions are connected via the crossing relations 



phJi 



-\v(v+u>)(v' —v) 



ll/[v+Cj)[v' — v) 



(B9a) 
(B9b) 



Therefore the only "new" (independent) quantity in the 
transverse (vertical) channel is (Eq. (|E>8b|) ) which 

corresponds to T ph xj- via the crossing relation Eq. (|B2|1 . 
Hence, in the following we will discuss only Eq. (|B8bj) in 
more detail: First of all we can perform the transforma- 
tion Wi = v\ — v of the summed index yielding 



x R 



n 



(BIO) 



In the next step we introduce the transformation v — > 
v. v' — > v + uj and uj — > v' — v and make use of the 
SU(2)-symmctry relation jD6| F$ +u){v '~ v) =-(F]g' u - 
P^' u ) = -F%? u . Furthermore we define f vv>ul = 



-,i/(i/+a;)(i/— u) 



Hence Eq. (|B10|1 reads 



puu'uj = {Wo, + ^^"G^JG^ + u)F£»'". 

(BU) 

This is exactly the same equation we derived for " 
(Eq. IB4b[) which means that 

f uv ' u = T%' u . (B12) 

Together with the definition of T this yields 

T™£" = _r^"+ w )("'-"). (B13) 

Hence, the transverse channel does not provide any 
"new" information (in the SU(2)-symmetric case), and 
T m and are, in fact, the only two independent func- 
tions for the two irreducible particle-hole channels. 



c. The particle-particle channel 

The particle-particle channel is completely indepen- 
dent of the two particle-hole channels and fulfills a cross- 
ing relation itself (Eq. (HH). The Bethe-Salpeter equa- 
tions for the three possible spin-combinations shown di- 
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FIG. 24. Bethe-Salpeter equations in the particle-particle 
channel. 



T^UU UJ F^U UJ I pUU UJ 

r t(riplet) ■— r pp,ti ~ l ~ pp,W 



(B15b) 



which correspond to Eqs. (fT7)l for the T's in complete 
analogy to the definition of the density and magnetic 
channel in Eqs. (|B3p . By adding and subtracting Eqs. 
(IB14bp and (|B14c[) one gets the Bethe-Salpeter equation 
for the singlet and the triplet channel 



F* 



= r u 



l i x - 

2^ 
1 1 



(B16a) 



F?'" = rr'" " 2 g E rr iiy '"G(^)G( W - Vl )F t 

(B16b) 

Writing the crossing relation (|D4b() in particle-particle 
notation yields 



t^vu'uj j-vv'iu—v—v') 



pp,tf 



jpv(u> — i'')(v'—v) _ j^v{uj — v')uj 
~ _ ph,tl ~ ~ W,tT 



pvvuj_ _ p> 



vv' {uj — v — 1>') _ jpv(uj — 1/ — v) _ t^v(u] — v')u 



(B17) 



Applying these relations to the definitions of singlet- and 
triplet-channel gives 



agrammatically in Fig. [M]read as 
1 1 



PP.Tt ~~ PP,tt 



rpw uj y\i/v uj 



(B14a) 



1 1 



X F 



v(uj — ui)uj 



PP,t4-,CTl(-<Tl) 



PP.Ti _ PP,T1 2/3^—^ PP,o-i(-cri).it 



x F 



v{uj — i^i)cj 



PP,t4.,0-l(-(Tl) 



(B14b) 
G(^i)G(w-i/i)x 

(B14c) 



The factor i appearing in these equations is needed in or- 
der to avoid double-counting since we are dealing with in- 
distinguishable particles (see e.g., [Uj]). The minus-sign 
in the reducible part again can be inferred from compar- 
ison with 2 nd -order perturbation theory. 
We see that in the particle-particle channel the ^t - 
vertex is completely independent from the two other spin- 
combinations, while T pp ^^ and T pp are not. Since they 
are coupled in the same way as T p h^ and T p h^i they can 
be decoupled introducing the linear combinations 



p^(cj — v' )uj pr/f i 



^y{uj~- v' )uj 



j^uv UJ 
1 t 



(B18) 



Inserting the crossing-relations for the ft - the singlet- 
and the triplet-vertex in Eqs. (|B14a[) and (|B16p yields 
again the standard matrix multiplication-form of the 
Bethe-Salpeter equations. Furthermore, combining these 
equations with the definition of the susceptibility in Eq. 
([9]) yields the corresponding Bethe-Salpeter equations for 
the generalized susceptibilities \ which read as 



uj -,vv ui 

As — AO 



1 1 



VV\U1 WiUNpl/rtU 1^1/ UJ 

2 ^2 i VAO, pp As / t AO, pp 



Xo,- 

(B19a) 



At — AO.pp 2 «2 VAO, pp 1" At ) L t A0,pp ) 

(B19b) 

where x s and \t ar< 3 defined analogously to the F , s in 

Eqs. (jBT5a|) and (|B15b|) . 

Solving Eqs. (|B19|) for T v / U and Y v t v ' u yields 

rr'" = - xo^p)- 1 + 2xoip]^'" (B20) 

r t — = /3 2 [4( Xt + xo,pp)' 1 -2 X -j p ]^. 



Considering the crossing relations ([B17)) and the SU(2) 
symmetry (Eq. ID6[1 one can express the singlet and the 



TpW UJ 

s(inglet) 



puis uj pvJJ UJ 

PPM pp,W 



(B15a) 
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triplet channel in the following way 



Appendix C: Parquet equations 



jpvv UJ J? vv UJ I <~>vv u 

T?VV UJ J? 1 ' 1 ' UJ 

r t - r PP ,n- 



(B21) 



This means that in the SU(2)-symmetric case there are 
only two independent irreducible particle-particle ver- 
tices namely, Y s and T t or Tfi and T^. 
However, this is to be expected since in the particle- 
particle case the and fi are connected via the crossing 
relation (|B17p . Because of that there is another possibil- 
ity to decouple the H from the ^-channel. Using the 
crossing relation (|B17|) we can eliminate T^^- from Eq. 

(|B14bp and obtain an equation containing only 



TpVV UJ 

pp,U 



(B22) 

Note that the factor ^ and the spin-summation have dis- 
appeared in this equation. Physically this result can be 
understood in the following way: The factor i was intro- 
duced in the particle-particle channel to avoid double- 
counting of diagrams since the two particles are indistin- 
guishable. This clearly holds for the ft -case - However, in 
the t|-case the spin can be fixed (i.e., no spin-summation 
in the Bcthc-Salpctcr equation) and hence, the two par- 
ticles are now distinguishable by their spin. 
Finally we write Eq. (|B22[) in terms of the corresponding 
susceptibility Xpp"i 



VV UJ 

x-pp,n 



v(uJ — V\)UJ y^V\V 2 UJ V 2 v' 



pp 



XppM 



ii Xo 



pp 



(B23) 

In contrast to Eqs. (|B19|) this equation does not have 
the form of a matrix-multiplication since it contains 

Xpp^ 1 ^ instead Xpp 1 ^ inside the sum. Nevertheless, 
it is possible to rewrite it by means of the substitution 
v' — >u— v' and the transformation v% — >uj— of the sum- 
mation variable v 2 . Considering that XcTpjT 2 ^" -1 ' ^ = 



In this section we give the explicit form of the parquet 
Eq. (|12p taking into their frequency dependence in terms 
of the density, magnetic, singlet and triplet channel in- 
troduced in the previous section. In order to simplify the 
notation we use the definition of reducible vertex $ 



fovv UJ T? vlJ UJ j^vv uj 



Hence, the parquet equations read 



d. m, s, t 



(CI) 



KVv'uj _ yw'uj , 3 ,„(„+„)(„'_„) 

lv d ~ L d +2 d 2 



}_^vv'(v+v'+uj) _ 3 ^vv'(v+v'+uj) 

2 s 2 4 



(C2) 



A vv'uj _ yvv'uj , ^faViu+u^iu' -u) _ }_tf ) v(v+uj)(v' -v) i 
m rn ' 2 9 



-|- (vW +w) _ Ijjjff'CH-f'+a;) 



a vv'uj _ yvv'uj _ }_*vv' (uj-v-v') 3 fa VV ' (tjJ-V—v') 



2~ d ' 2 

}_Jf ) v {'^- v '){ v '- v ) , ^_^v(uj-v')(v' -v) 

2 d 2 " 



(C4) 



j^vv'uj -pvv'uj _!_ cE) 2 ^' u ~ v ' ^ ^_fi*vv' (uj — v — v') 

t t 2 d 2 



(C5) 

For the A s and A t particle-particle notation was adopted. 
Since at the level of A no dependency on an irreducible 
channel (ph, ph or pp) is present A s and At can be ex- 
pressed in terms of the A^ and A rn 



A VV'UJ _ 1 A VV(UJ — V — V') _ 3 . V v' (uj-V—v') 
~ r,^d 2 m 



A'; 



1 . vv(uj-v-v') 1 . vv'(uj-v-v') 



(C6) 



v(uj — v')u 

Xpp^l 



V vvxuj v(uj — vi)u\ i j^v 1 {uj — V2)uj^v 2 v'uj 

^2 yXo.pp — Xpp,tl 



J )r 



PPM 



r>p 



With the definition Xpp "i = Xpp^i' (and the same for 



(B24) 



the T's) one gets the Bethe-Salpeter equation (|B24|) in 
the usual form of a matrix multiplication 



-w'uj J_ \^(-,V 

Xpp^i ^2 / AXp 



VV\UJ 

pp 



-vvxuj \-pVtv 2 uj v 2 v uj fnoi:\ 
XppAl' PPM Xo,pp i 



It can be solved for f yielding 

T ppM = P 2 [(XppM ~ Xo.pp)" 1 + Xo 



X I VV UJ 

pp\ 



(B26) 



Appendix D: Symmetries 

In this appendix, we summarize for convenience some 
symmetry properties of one-particle Green's function G 
and the generalized susceptibility x- 



1. Time- Reversal Symmetry 

A system without spin-orbit coupling is invariant un- 
der time-reversal if its Hamiltonian T-L (assumed to be 
time- independent) is a real function of the momentum 
operator p and the position operator x. This usually 
holds in absence of an external magnetic field. It can be 
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shown that one can always find real eigenfunctions ip(r) 
in this case and and analogously, the n-particle Green's 
function is a purely real function of the (imaginary) times 



These relations can be easily proven by rotating all spins 
through an angle tt about the x- or y-axis. Furthermore, 
performing a rotation through an angle i.e., rotating 
a spin in z-direction into the a;y-plane, yields 



Gn(ri, . . . ,r 2n ) = G„(ri, . . . ,T% n ). 



(Dl) 



This property of the n-particle Green's function can be 
easily proven by passing on to its functional integral 
representation^. 

As the AIM defined in Eq. ([2]) complies with all above- 
mentioned conditions, the imaginary-times 71-particle 
Green's functions are real. Hence, one can derive the 
following relations for the one- and the two-particle 
Green's functions (i.e., the generalized susceptibility) of 
this model in frequency-space 



G» = G{-v) 



Xc 



Act'ct ' 



(D2a) 



(D2b) 



Let us also give an equation relating the generalized sus- 
ceptibility x to its complex conjugate 



(x 



<7<7' ) 



{- v ')(-v){-u) 



(_„) (-</)(-o;) 



(D3) 



2. Crossing Symmetry 

This symmetry is simply a consequence of the Pauli- 
principlc, i.e., exchanging two identical fermions leads 
to a minus-sign in the wave function. Considering Eq. 
(|7ap the exchange of annihilation operators in the time- 
ordered matrix element yields a minus-sign and leads to 
an exchange of the corresponding frequencies v' and v+lo. 
Taking into account additional xo-contributions one gets 
the following crossing relations for \, F and A in particle- 
hole notation 



^.Vv'u _ r v(v+u)(v' -v) 



v(v+w)(v'-v) uv'ui 
X-oo' ' AO i 



crcr' 



crcr' 



_pv(y+u)(y' —v) 



_^v(v+u})(v' — v) 



(D4a) 
(D4b) 



(D4c) 



3. SU(2) Symmetry 

If the Hamiltonian of the system does not contain 
terms breaking rotation symmetry (e.g., a magnetic 
field), the x's and the F's satisfy some specific relations. 
Every matrix-element hast to fulfill spin-conservation, 
e.g., Cr-M = 0. The one-particle Green's function is in- 
dependent of the spin, i.e., = G44 = G. At the two 
particle level similarly xtf = Xu and Xn = Xlt nold - 



VV u 



vv'u) _ v (v+u)(v'-v) 
Xa(-a) Affj-ir) 



Xo 



^acr - t <r(-<r) 



F 



cr(-cr) 



(D5) 
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4. Mapping onto the Attractive model 

Usually partial particle-hole transformations which 
map repulsive onto attractive Hubbard interactions are 
defined for lattice systems^. Obviously, one can find an 
equivalent (but local) transformation for the correspond- 
ing AIM as will be shown in the following. 
Starting point is the Hamiltonian of the AIM (Eq. ((5J)) 
containing also the chemical potential term — /x(n-|- + 
since we consider a grand canonical ensemble 



(D7) 



The sum over £ (bath sites) ranges from 2 to N, £ = 1 

denotes the impurity (i.e., a^} = c# ). The (partial) 
particle-hole transformation we are considering is defined 
by the unitary operator W 



(D8) 



The action of the transformation W on the creation and 
annihilation operators is given by 



W\c\,ci)W = (-1)^(01,01) 
VVt(c|,c t )VV = (c|,6 t ) 



VV f (ai 



VVt(a] a it )W = (-l) N (al f ,a it ). 



(D9a) 



(D9b) 



This means that for a =\. the annihilation and creation 
operators are interchanged, while the t-°P era tors are 
not modified by the transformation VV (despite a phase- 
factor (— 1) N ). Therefore VV is coined partial particle- 
hole transformation. 

The transformation of the AIM- Hamiltonian given in Eq. 
(lD7l) yields 



IN 



W 



£=2 



(D10) 
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We are now restricting ourselves to an even number of 
bath sides, i.e., N has to be odd. Furthermore, we as- 
sume that the bath levels are distributed symmetrically 
around 0, i.e., eg = — s e+ N for I = 2..^ + 1. In addi- 
tion, the hybridization between the bath and the impu- 
rity should be the same for positive and the correspond- 
ing negative bath-energies which means that Vt = V i+ N_ 

for I = 2.. y + Hence, the negative energy-sector of the 
bath is completely equivalent to the positive one. Per- 
forming the index-transformation I <-> [l + ^) for the 
4-spins in Eq. (|D10[) changes the minus-sign in front of 
ega\^aei back into a plus-sign as in the original Hamilto- 
nian. 

Furthermore, choosing the chemical potential as 

^ = §, (DH) 

in Eq. (|D10|) one retrieves the same structure of the 
original Hamiltonian whereas only the sign of U has 
changed (the constant contribution —fj, + ^^=2 £ ? can 
be neglected). 

This way it has been shown how the transformation W 
maps the repulsive AIM Hamiltonian (£/>0,/i = -j)on 
the attractive one (— U, /x = — ~) provided that the ad- 
ditional conditions for N, si and Vi are fulfilled, which is 
the case for the particle-hole symmetric AIM associated 
to the DMFT solution of the half-filled Hubbard model 
considered here. 

Next, we discuss some symmetry-relations for the n- 
particle Green's functions G^ a CT2ti and the general- 
ized susceptibility xWaa'- The additional index U indi- 
cates whether the quantity under consideration is calcu- 
lated for repulsive (U or +U) or for the corresponding 
attractive (-U) model. 

First, applying the particle- hole transformation to the 
one-particle Green's function with spin-f lets this func- 
tion unchanged, since W only acts on the ^.-creation- and 
annihilation operators, i.e., 

G^(t 1 ,t 2 ) = G\- u \t 1 ,t 2 ). (D12) 

For the Green's function with spin-J, 64, and ct change 
their role which leads to an exchange of T\ and t 2 as well 
as to an additional minus-sign 

G u u {t 1 ,t 2 ) = -G{- U \t 2 ,t 1 ). (D13) 

For the SU(2) symmetric case = G44 = G one can 
combine relations (|D12[) and (|D13jl and gets 

G u (n,T 2 ) = -G u (T 2 , n ). (D14) 

which means in Fourier-space 

G*{v) = -G{v), (D15) 



expressing the fact that in the particlc-holc symmetric 
case the one-particle Green's function is purely imagi- 
nary. 

Taking the limit t 2 — > ti+ (i.e. t 2 — > ti and t 2 > t±) in 
Eq. (|D14[) leads to the result that the average density at 
the impurity (n) = n = 1, which means that the system 
is " half-filled " in the particle- hole symmetric case. 
Next, we consider the two-particle Green's function, i.e., 
the generalized susceptibility. As in the onc-particle-case 

G^ tttt (Ti,T2,T3,T 4 ) = G 2t ^.(n,T2,T 3 ,T4,). (D16) 

The two-particle Green's function containing only .[-spins 
transforms under VV as follows: 

G 2,UU ( r i , r 2 , r 3 , r 4 ) = Gjj^+U ( T 4 > T 3 > T 2 , n ) . (D 1 7) 

Combining cqs. (|D16p and (|D17|) and using again SU(2)- 
symmctry yields 

G 2j ^^.(ti,T2,T3,T4,) = G 2j -fftf.(T4,,T3,T2,Ti). (D18) 

In Fourier-space this relation states that the two-particle 
Green functions are purely real and the same holds true 
also for the susceptibilities 

Furthermore we want to study how the ^-function trans- 
forms under the particle-hole-transformation. In the 
corresponding matrix element only the operators corre- 
sponding to the times T3 and T4 carry ^-spins and there- 
fore VV acts only on them 

In Fourier-space this is equivalent to the transformations 
(z/ + oj) -> (— v') and v' -> (— v' - w), i.e., 

VV ' U> V{v+Uj){v' — V) v{ — v' ~Uj)uj v{v-\-Ld){ — V — V — Oj) 

Xu,n~*o — — X(_j/),-f4. +*o 

(D21) 

Using SU(2) symmetry on the left hand side of this equa- 
tions yields 

_ v(v+u)(v' -v) _ v{-v'-ui)u v{v+u){-v-v' -u) 
X-U,m — X(-U),1l ^ Xo 

(D22) 

Performing the frequency transformation v' — > v — u> and 
uj — y v 1 — v and transforming the right hand side to the 
particle-particle notation gives 

vv'(-ui) _ v(u-v')u> vv'u IV\OX\ 

X.U,m ~ X-(-U),pp,ti. Xo, P p ■ yuto) 

This equation can be interpreted as follows: The inver- 
sion Xrjm ^ yields Tm as discussed in Sec. |5] 
The inversion of the quantity on the right hand side of 
Eq. (|D23[) gives the irreducible fi-vertex in the particle- 
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particle channel, i.e., ^\-u)wn ( see Eq ' E2S|. Hence, which is also shown diagrammatically in Sec. HVl 

If one performs the sum over v and v' in Eq. (|D23|) one 
piV(-w) _ j,v{u)-u')u) (D24) sees fluctuations of the spin for the repulsive model 

m (-u),pp,W are mapped on fluctuations of an electron-pair for the at- 

tractive case. This is consistent with the well-known fact 
that for a lattice model the anti-ferromagnetic instability 
for U > corresponds to the superconducting instability 
in the attractive model. 
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